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Optomechanical devices sensitively transduce and actuate motion of nanomechanical structures 
using light. Single-crystal diamond promises to improve the performance of optomechanical devices, 
while also providing opportunities to interface nanomechanics with diamond color center spins and 
related quantum technologies. Here we demonstrate dissipative waveguide-optomechanical coupling 
exceeding 35 GHz/nm to diamond nanobeams supporting both optical waveguide modes and me¬ 
chanical resonances, and use this optomechanical coupling to measure nanobeam displacement with 
a sensitivity of 9.5 fm/\/Hz and optical bandwidth > 150nm. The nanobeams are fabricated from 
bulk optical grade single-crystal diamond using a scalable undercut etching process, and support 
mechanical resonances with quality factor 2.5 x 10® at room temperature, and 7.2 x 10® in cryogenic 
conditions (5K). Mechanical self-oscillations, resulting from interplay between photothermal and 
optomechanical effects, are observed with amplitude exceeding 200 nm for sub-/rW absorbed opti¬ 
cal power, demonstrating the potential for optomechanical excitation and manipulation of diamond 
nanomechanical structures. 


Nanophotonic optomechanical devices w allow 
chip-based optical control of nanomechanical resonances 
with a precision reaching the standard quantum limit lu¬ 
ll] , enabling tests of quantum nanomechanics [7] I^HTT] , as 
well as technologies for sensing □[EHIII and information 
processing [ISHIT]. Single crystal diamond devices are 
particularly exciting for these applications. In addition 
to possessing desirable mechanical and optical properties, 
diamond hosts color centers [18] whose highly coherent 
electronic and nuclear spins are promising for quantum 
technologies. Recently, piezoelectric-actuated nanome¬ 
chanical manipulation of electronic spins in diamond 
color centers [I9U2T| and quantum dots O [23] have illu¬ 
minated the possibility of using diamond nanomechanical 
devices |71IE5] to mediate phonon-spin interactions at 
the quantum level. Adding optomechanical functionality 
to these devices will allow efficient coherent optical ma¬ 
nipulation and readout of the nanomechanical stuctures, 
and provide a path towards implementing an efficient in¬ 
terface between optical, mechanical and electronic quan¬ 
tum systems. Such an interface would play a crucial role 
in many proposed studies, including normal mode cool¬ 
ing |29l 130] , quantum nonlinear optomechanics m, and 
spin-squeezing [32]. 

Here we demonstrate a diamond waveguide optome¬ 
chanical system providing sensitive and broadband read¬ 
out of ultrahigh mechanical quality factor (Qm) reso¬ 
nances which can be excited into large amplitude optome¬ 
chanical self-oscillations, revealing the nonlinear nanome¬ 
chanical properties of the nanobeams. Dissipative op¬ 
tomechanical coupling in this system relies on efficient 
phase-matched evanescent coupling between a hber ta¬ 
per waveguide and optical modes of ultrahigh-Q^ dia¬ 
mond nanobeam resonators. The nanobeams are fab- 
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ricated from single-crystal diamond chips using oxygen 
plasma undercut etching. They support mechanical res¬ 
onances whose motion is transduced by the optical fiber 
taper-nanobeam coupling, resulting in an optomechan¬ 
ical interaction characterized by an effective dissipative 
optomechanical coupling coefficient of 35 GHz/nm, and 
a displacement sensitivity of 9.5 fm/VHz. In combina¬ 
tion with low nanobeam mechanical dissipation {Qm > 
7.2 X 10®), interplay between pN photothermal forces and 
optomechanical coupling is shown to excite nanomechan¬ 
ical self-oscillations with amplitude > 200 nm, which 
renormalize the nanobeam mechanical frequency and 
provide a measure of the nanobeam’s internal stress and 
buckling amplitude. The excellent displacement sensi¬ 
tivity, broadband coupling, and ultrahigh-Q^ demon¬ 
strated here will enable diamond nanobeam waveguide 
optomechanical implementations of sensors [injus] and 
quantum nanomechanical systems 011 ], while the ability 
to optomechanically excite nanomechanical resonances 
shows the potential for optical manipulation of local 
stress fields and accompanying phonon-spin coupling in 
diamond color centers [TMHI33H35], and realization of 
hybrid photon-spin-phonon systems m- 


I. WAVEGUIDE OPTOMECHANICS 

Experiments in diamond nanophotonics and nanome¬ 
chanics have recently been advanced by the availability of 
high-quality diamond chips grown using chemical vapor 
deposition. While optomechanical devices can be fabri¬ 
cated from polycrystalline diamond films m. single crys¬ 
tal diamond thin hlms which retain the desirable com¬ 
bination of optical, mechanical and quantum electronic 
properties introduced above are not commercially avail¬ 
able, and must be manufactured using ion-implantation 
[55H^ or wafer-bonding [^[2311331110] techniques. Fab- 
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FIG. 1. Diamond nanobeam waveguide optomechanical system, (a), Scanning electron micrograph (SEM) images of a single 
crystal diamond nanobeam-waveguide. Dark high-contrast regions are due to variations in thickness of titanium deposited for 
imaging purposes, (b), Schematic of fiber taper evanescent waveguide coupling geometry and illustration of the waveguide- 
optomechanical coupling process. Nanobeam mechanical resonances, whose typical displacement profiles are shown, modulate 
the distance between the input waveguide, changing K{h). (c). Dispersion of riefj of the waveguide modes of the fiber taper 
(diameter 1.1 /rm) and diamond nanobeam {wxd — 460x 250nm^) when they are uncoupled (dashed lines, red: fiber taper, blue: 
nanobeam) and coupled with h — 300 nm (solid lines). Shaded region indicates the laser scan range used in the experiments. 
Insets: mode profiles {x — y plane, dominant electric field component \Ey\) of the symmetric and antisymmetric coupled modes. 


ideation of devices directly from bulk single crystal dia¬ 
mond chips is desirable, and until now has relied on either 
ion beam milling mi: or less damaging and more scal¬ 
able inductively coupled plasma reactive ion (ICPRIE) 
Faraday cage angled-etching [281 Il2]- In this work we 
demonstrate an ICPRIE process which does not require 
a Faraday cage, and instead utilizes diamond undercut 
etching to fabricate nanobeams, such as those shown in 
Fig. 0 a), from bulk single crystal diamond. This process, 
summarized in Fig. and discussed in more detail in Ap¬ 
pendix is inspired by earlier bulk silicon nanofabrica¬ 
tion |43| and relies on an oxygen based ICPRIE process 
operated with zero-RF power, high-ICP power and an 
elevated sample temperature (250 °C) to etch diamond 
quasi-isotropically along diamond crystal planes. The 
process is uniform over the diamond chip area and is fully 
compatible with standard commercial etching tools. 

Using this diamond undercutting process, nanobeams 
were fabricated with widths w = 300—540 nm and 750 nm 
and lengths L = 50 — 80 /im on the same chip. Nanobeam 
thickness d was adjusted by controlling the etching times 
together with the nanobeam width w. For the etch pa¬ 
rameters used here, the nanobeam thickness was d ~ 
250 — 350 nm, depending on w. Owing to the crystal- 
plane sensitive nature of the undercut etch, for the etch 
parameters used here the narrow nanobeams have a flat 
bottom surface, while the w = 750 nm nanobeams have 
a triangular bottom surface mimicking the ridge visible 
in Fig. Ba). These bottom surfaces can be flattened by 
longer undercut etching. 


The fabricated nanobeams support mechanical reso¬ 
nances, the lowest order of which are illustrated in Fig. 
Bb), with frequencies and effective mass in the MHz 
and pg range, respectively. The nanobeams also sup¬ 
port waveguide modes, which can be evanescently cou¬ 
pled with high efficiency to an optical fiber taper [44j . 
as shown schematically in Fig. [^b), by tuning the de¬ 
vice geometry to match the nanobeam mode phase ve¬ 
locity with that of a fiber taper mode. Despite the re¬ 
fractive index difference of the Si 02 fiber taper and di¬ 
amond, phase-matching is realized in nanobeams with 
subwavelength cross-section. This is illustrated in Fig. 
j^c) which shows the effective refractive index dispersion, 
neff(A), of the fundamental TE modes of a fiber taper and 
of a diamond nanobeam, together with that of the “su¬ 
permodes” of the evanescently coupled waveguides (n^) 
[45] . Phase-matching occurs at Ao ~ 1570nm, where 
an anti-crossing in n^(A) indicates that the waveguide 
modes are coupled. 

Evanescent waveguide coupling, which has been well 
studied in the context of photonic and optoelectronic 
devices |45l 146] . is shown here to provide sensitive op¬ 
tomechanical readout of mechanical fluctuations which 
modulate the waveguide separation, h. In contrast to 
cavity-optomechanical systems which exploit dispersive 
coupling between mechanical fluctuations and narrow- 
band optical resonances m, waveguide optomechanical 
coupling utilizes wide-band dissipative optical transduc¬ 
tion of the nanobeam mechanical position. The sensitiv¬ 
ity of the evanescent waveguide optomechanical coupling 
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FIG. 2. Process flow for creating diamond nanobeams using quasi-isotropic reactive ion undercut etching. 
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FIG. 3. Efficient evanescent waveguide coupling, (a), Fiber 
taper transmission T{X) vs. wavelength for h = 200 nm. In¬ 
set: variation of T with nanobeam width w and A for approx¬ 
imately constant h. (b), Dependence of fiber taper trans¬ 
mission on fiber height. Experimentally measured values for 
To (red points) are fit with the model described by Eq. Q 
(red line). The corresponding K{h) extracted from the fit 
(blue line) and predicted from simulation (blue points) are 
also shown. 


can be derived from coupled-mode theory (Appendix [b|), 
and is determined by the dependence of the normalized 
fiber taper transmission, T on the separation h between 
the waveguides: 

T{h,X)=cos^{sL,)+(^^^ sm^^^ 

Here L^. is the coupler interaction length determined by 
the fiber taper geometry, A/3(A) = Aneff(A)27r/A is the 


propagation constant mismatch of the uncoupled waveg¬ 
uide modes, + Ap'^jA. K{h) is the per-unit 

length amplitude coupling coefficient, which depends on 
the overlap of the evanescent fields of the coupled waveg¬ 
uides. The theoretical displacement sensitivity achiev¬ 
able by monitoring fluctuations in T is 


s:if) = 


2T{h,X)PdfuJo/ri' 


qe 




Pj (dT/dn) (dhi/dh) (dh/dx) 


( 2 ) 


where is the single-sided noise equivalent optical 
power spectral density (units of W^/Hz) of the detec¬ 
tor and technical noise, and Pd is the detected out¬ 
put power in absence of coupling (T = 1). The first 
term in the numerator accounts for photon shot noise, 
where is the detector quantum efficiency. This ex¬ 
pression neglects radiation pressure backaction, which 
is not significant for the measurements presented here, 
but will ultimately limit s°. The impact of coupler 
geometry, operating condition, and waveguide design 
on detection sensitivity is described by the denomina¬ 
tor of Eq. ([^. Sensitivity is maximized for mechani¬ 
cal resonances whose displacement x efficiently modu¬ 
lates h, i.e., \dh/dx\ ~ 1. Strong evanescent overlap 
enhances |9 k/9/i|, while phase-matching and operation 
near K{h)Lc ~ tt/A,3tt/A, ... maximizes |9T/9 k|. 

Experimental observation of the optomechanical prop¬ 
erties of the diamond waveguide optomechanical system 
was performed by measuring T{t) = T + ST(t) of a dim¬ 
pled optical fiber taper [33] positioned in the nanobeam 
near-held for varying h and A. The efficiency of the 
evanescent coupling determines the time-averaged trans¬ 
mission, T. Eluctuations of the nanobeam position, to- 

S ither with other noise, are imprinted on ST{t). Eigure 
a) shows T(A) when the flber taper is positioned at 
h ~ 200 nm above the center of a nanobeam (w x d = 
460 X 250 nm^). Minimum transmission To = 0.05, corre¬ 
sponding to coupling efficiency of 1 — To = 95% assuming 
negligible Insertion loss, is observed near Ao = 1560 nm. 
The 3 dB bandwidth AA > 150 nm is consistent with pre¬ 
dictions from the n^(A) anti-crossing in Fig. He), and 
increasing w was observed to increase Ao (Fig. p a) inset). 
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(b) Frequency [MHz] (d) 



FIG. 4. Nanomechanical and optomechanical properties, (a), Measnred Sx{f) when the fiber taper is conpled to a nanobeam 
{L X w X d = 60 X 0.46 x 0.25 /rm®) in ambient conditions, with Pd ~ 100 fj,W. The vertical axis and noise-floor (dashed line) 
of Sx = 9.5 fm/\/Hz were calibrated from the fit to the / = 1.3 MHz mechanical resonance thermomechanical displacement 
spectrum (red solid line), (b), Top: observed (points) and fit with a function oc \dT/dh\ (solid line). Bottom: 

displacement sensitivity s'^{h) of the fundamental out-of-plane resonance from measurement (blue points) and predicted from 
the measured T{h) for varying input power (solid lines). The starred point corresponds the sensitivity and operating condition 
of the measurement in (a). Measured Sv{f) for a nanobeam with Lxwxd — 60 x 0.75 x 0.3 /rm® (c), in vacuum (d), at 5K. 
Insets: fits to observed mechanical resonances, and corresponding Qm- 


consistent with expected behavior. 

The coherent nature of the waveguide coupling is re¬ 
vealed by To{h). As shown in Fig. ^b), To{h) is mini¬ 
mized at ~ 200 nm, where kL^ = 7r/2. For h < 200 nm. 
To increases with decreasing h due to the co-directional 
coupling undergoing more than a half “flop” and light 
coupled into the nanobeam being out-coupled back into 
the fiber taper [45]. In contrast, incoherent scattering 
loss increases monotonically with decreasing h |48j . and 
is small in the system studied here. As shown in Fig. 
ib), Eq. (1^ fits To{h) well with an exponentially de¬ 
caying K{h)as a fitting parameter, which in turn agrees 
closely with k(/i) predicted from coupled-mode theory 
for interaction length Lc ~ 7 /im. This Lc is consistent 
with the observed curvature in optical images of the fiber 
taper dimple. 

The optomechanical properties of the coupled waveg¬ 
uides were observed from the power spectral density 
Sv{f) of the photodetected signal generated by fluc¬ 
tuations in output power Pd5T{t). Figure |^a) shows 
the equivalent displacement spectral density Sx{f) (units 
m/VHz) of the nanobeam motion when the fiber taper is 
positioned h ^ 100 nm above the nanobeam. Peaks from 


thermally driven nanobeam resonances are clearly visi¬ 
ble at frequencies /„ = [1.3, 3.1,4.2] MHz, correspond¬ 
ing to the fundamental out-of-plane (vi), fundamental 
in-plane (hi), and second order out-of-plane (V 2 ) reso¬ 
nances, whose simulated displacement profiles are shown 
inFig.Qb). Resonance labels can be determined by com¬ 
parison of fm with simulations and by measuring trans¬ 
duction sensitivity as a function of in-plane taper posi¬ 
tion. Sx{f) was obtained by calibrating Sv{f) to the the¬ 
oretical thermomechanical power spectral density of the 
Vi resonance (m = 7.6 pg), as described in Appendix [Pj 
The observed technical noise floor of s° = 9.5fm/Hz^' 
is below that of other broadband integrated waveguide 
measurements ETlIlHlIlHj, and is more than an order of 
magnitude more sensitive than typical free-space reflec¬ 
tion measurement techniques |28j . This sensitivity can 
be further improved to the sub-fm/Hz^^^ range by in¬ 
creasing Lc and Pd- 

The h dependence of and s° for the Vi res¬ 

onance, shown in Figure |^b), provides a direct measure 
of the optomechanical coupling. The on-resonance signal 
is directly related to the slope of T{h): Sy{h; oc 

\dT/dh\, with a distinct minimum when kLc = 7r/2. 
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Highest sensitivity is observed at h = 100 nm {kLc = 
37 r/ 4 ), and s°{h) agrees well with predictions from the 
experimentally characterized T{h) and Eq. (H- The min¬ 
imum s%{h) in Fig. |^b) is degraded compared to the 
sensitivity of the measurement in Fig. |^a) due to op¬ 
erating at lower Pd- dT/dh is related to the effective 
dissipative optomechanical coupling coefficient of the sys¬ 
tem, defined by the change in waveguide tunneling rate 
of photons in the coupling region for a given change in h: 
ge j27r c |^|/(27rngLc-\/?o) > 35GHz/nm where c/rig 
is the group velocity of light in the coupler. 

Nanobeam mechanical dissipation in the measure¬ 
ments discussed above was dominated by damping from 
the ambient air environment. To assess the nanobeam 
mechanical properties, measurements were performed in 
vacuum and low-temperature conditions. Generally, af¬ 
ter reducing air damping was observed to increase 
above 10 "* in all device geometries, with the lowest dissi¬ 
pation {Qm > 10 ®) observed in the larger cross-section 
(w X d = 750 X 300nm^) devices. Figures l^c) and 
®d) show mechanical resonances of an ultrahigh-Qm 
nanobeam measured in vacuum (80 ^Torr) and in cryo¬ 
genic conditions (T = 5 K and 5 /rTorr). In vacuum, the 
vi and hi resonances have Qm ~ 1-4 x 10® and 2.0 x 10® 
respectively. At low temperatures, dissipation was fur¬ 
ther reduced, such that Qm ~ 1-7 x 10® and 7.2 x 10® 
for the Vi and hi modes, respectively. These values of 
Qm are higher than previous reports of nanomechanical 
devices fabricated from single crystal optical grade dia¬ 
mond [251 IMl HH] ■ All measurements were made at low 
power to avoid inducing optomechanical linewidth nar¬ 
rowing. Fabricating devices from electronic grade ma¬ 
terial or with larger dimensions may allow for further 
increased in Qm, as in Ref. [2H] . 

Given the demonstrated device performance, the mea¬ 
surement precision required to reach the standard quan¬ 
tum limit is = \/hQm/{muJm) = 3.1 fm/vTdz where 
m = 23 pg is the effective mass of the hi resonance 
[5]. While this is ^ 3 times smaller than the mea¬ 
surement sensitivity demonstrated here, promising ap¬ 
proaches for reducing s° below include improving 

the optomechanical coupling by increasing Lc, and fabri¬ 
cating higher-Qm devices from electronic grade diamond 
material [26]. 


II. TUNABLE NONLINEAR DYNAMICS 

Backaction from waveguide optomechanical coupling 
can dramatically modify the nanobeam dynamics, and is 
shown here to amplify the nanobeam motion and reveal 
nonlinear nanomechanical properties of the device. These 
effects can be tuned by adjusting the waveguide position, 
and depend critically on the presence of internal stress in 
the nanobeam. Demonstration of optomechanically mod- 
ihed nanobeam dynamics is shown in Fig. [^a), where 
time-resolved T{t) of a fiber taper coupled to a high- 
aspect ratio nanobeam {L x w x d = 80 x 0.48 x 0.25 /rm®) 


in a vacuum environment is recorded while h is discretely 
stepped (30 nm/step, 0.8 steps/s. Pi ^ 300^W). For 
large h, T(h) behaves similarly to the ambient condition 
measurements in Fig. |^a). However, at h ^ 400 nm fluc¬ 
tuations 5T{t) become large-amplitude self-oscillations 
with a peak-to-peak change AT > 0.3, corresponding to 
nanobeam displacement exceeding 200 nm, as illustrated 
in the inset to Fig. [^a). To the best of our knowledge, 
this amplitude is larger than in previous reports of on- 
chip externally driven single crystal |50j and polycrys¬ 
talline |27j diamond nanobeams. Finite elements simu¬ 
lations predict a variation of axial stress at the center of 
the nanobeam on the order of 70 MPa for the vi reso¬ 
nance of the compressed nanobeam oscillating with this 
amplitude. When h < 225 nm, the self-oscillations stop. 
Figure [^b) shows a spectrograph of this data, where the 
vi resonance near fm = 430 kHz with Qm = 2.5 x 10^ 
is observed to increase in amplitude, shift to lower fre¬ 
quency, and generate harmonics. 

In contrast to diamond nanobeam oscillations driven 
by external actuation |50j or resonant optical modula¬ 
tion the observed self-oscillations are driven by 

the dynamic interaction between a delayed photother- 
mal force [HI US] and the waveguide optomechanical cou¬ 
pling. The resulting optomechanical backaction renor¬ 
malizes the mechanical dissipation rate from 7 ^ to 


Ir, 

Ir, 


f Qm 


T WmT dT 
k 1 +W^r 2 dx 


(3) 


where t is photothermal force response time, and P 
is the photothermal force strength per unit absorbed 
power Tabs = Q — T)Pi(Li, defined by the corresponding 
nanobeam deflection and spring constant k = muJm of the 
mode of interest. Li is the distance over which input light 
propagates in the nanobeam before being out-coupled, 
and C is the waveguide absorption coefficient per unit 
length. Figure [^c) compares experimentally observed 
'y'm{h) with theoretical values obtained from Eq. (§ in¬ 
put with measured T{h) as well as parameters discussed 
below. Predicted 7 m(h) agrees well with experiment, par¬ 
ticularly in reproducing the range of h over which the op¬ 
tomechanical coupling {\dT/dx\) is sufficiently large such 
for 7 m < 0, resulting in self-oscillations. Note that dis¬ 
agreement for h < 100 nm is in part related coupling be¬ 
tween fluctuating fiber position and Wm, which becomes 
significant compared to 7 m for small h (see Appendix [f|) 
In addition to exhibiting ultrahigh-Qm and strong dis¬ 
sipative optomechanical coupling, the nanobeams have 
two properties which make them sensitive to optome¬ 
chanical photothermal actuation. The nanobeam ther¬ 
mal time constant r ~ 0.7 /iS, calculated using finite el¬ 
ement simulations, is on the same timescale as Wm- As 
seen from Eq. (§, this is a necessary condition for sig¬ 
nificant feedback from the photothermal force [5IHS1|. A 
more subtle but equally important property is the pres¬ 
ence of compressive stress and accompanying buckling in 
the nanobeam, which as discussed below, dramatically 
enhances P. 
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FIG. 5. Nanobeam self-oscillations, (a), Time resolved T{h{t)) of a phase-matched fiber taper-nanobeam {L x w x d = 
80 X 0.50 X 0.25/rm®) system as h is reduced in discrete steps (visible as sharp steps in T). Inset: large amplitude oscillation of 
the nanobeam position, (b), Spectrograph showing the power spectral density of the data in (a), (c), Predicted renormalized 
mechanical dissipation rate compared with measured values, (d), Scatter plot of shift in nanobeam resonance frequency vs. 
oscillation amplitude squared for data in the inset regions from (b) at the onset of self-oscillation. The solid line is a linear 
guide to the eye. 


For the nanobeam geometry studied here, hnite ele¬ 
ment analysis presented in Appendix indicates that 
can be enhanced by over two orders of magnitude in a 
compressed and buckled nanobeam. The level of com¬ 
pressive stress in a nanobeam can be estimated from 
the deviation of ujm from the nominal value expected 
for an uncompressed device [28j . The nanobeam stud¬ 
ied in Fig. is observed to be in a post-buckled state 
with measured a;m/27r ~ 430 Hz significantly lower than 
the nominal value expected (680 kHz) for an uncom¬ 
pressed nanobeam. Matching hnite element simulated 
and observed uim predicts an internal compressive stress 
of ~ 37 MPa and an accompanying buckling amplitude of 
X = —122 nm. For this value of internal stress and buck¬ 
ling, a photothermal force of = — 26pN/p,W is pre¬ 
dicted, which is over 100 times larger than predicted for 
an uncompressed device. Note that the negative signs of 
T and x indicate that the photothermal force and buck¬ 
ling amplitude are in the down direction, consistent with 
all of the self-oscillation behavior discussed in this sec¬ 
tion. 

Given the above device parameters, the only free htting 
parameter needed to match Eq. ([^ with experiment is 


the waveguide absorption coefficient, which is found to be 
C = 0.12 cm“^, corresponding to an optical loss rate effec¬ 
tive quality factor of Qo 6.6x 10^. This absorption rate 
is consistent with loss in other diamond nanophotonic de¬ 
vices [42], and indicates that only a fraction (^ 100 nW) 
of the input power is absorbed and responsible for driving 
the self-oscillations. 


Nonlinear coupling between nanobeam oscillation am¬ 
plitude and uJm provides an additional probe of the in¬ 
ternal stress and buckled state of the device. Softening 
of the mechanical resonance frequency by AijJm{h) is ob¬ 
served at the onset of self-oscillations, as highlighted in 
Fig. ^h). Nanobeam softening and hardening is a well 
known indicator of internal stress [2711501, and the soft¬ 
ening observed here is found to be closely related to the 
buckled nanobeam state. Nanobeam buckling breaks the 
device vertical symmetry and introduces a nonlinear soft¬ 
ening term to the nanobeam dynamics which counteracts 
the hardening effect described by the intrinsic Duffing 
nonlinearity [55j . This competition between nonlinear 
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effects is given by, 

. , 3 _n 15 o, , , 

= —(otts - X (4) 

UJrn o 

where v is the oscillation amplitude, and is the Duff¬ 
ing coefficient of the unbuckled nanobeam as described 
in Appendix]^ Equation Q clearly shows the softening 
influence of x. The optomechanical system studied here 
provides a direct measurement of Awm and v for varying 
/i, allowing x to be estimated experimentally. Figurej^d) 
shows a scatter plot of measured Awm as a function of v'^ 
for varying fiber taper position. As predicted from Eq. 
Q, Aujm and are found to be linearly related, and 
from this data and Eq. Q, x = —98nm is estimated, 
which is in excellent agreement with finite element sim¬ 
ulation predictions of x given above. 

During the self-oscillation limit-cycle, a low funda¬ 
mental self-oscillation frequency of ^ 180 kHz is mea¬ 
sured, consistent with the behavior of a nanomechanical 
resonator oscillating between buckled states, as observed 
by Bagheri et al. |1^ . Further evidence of this behav¬ 
ior is found in Fig. ^a), which shows that during self¬ 
oscillations T decreases, indicating that the nanobeam, 
initially in a buckled down state, moves on-average closer 
to the fiber taper while self-oscillating. Note that har¬ 
monics which emerge during self-oscillation, becoming 
stronger and more numerous as v increases, result from 
both the nonlinear response of T{h) and nonlinearities of 
the nanobeam [56| . 


III. DISCUSSION AND CONCLUSION 

The waveguide optomechanical interface demonstrated 
here possess a unique combination of high sensitivity, 
broad bandwidth, high quality single-crystal diamond 
material, and high-Qm, and has potential to allow mea¬ 
surement of quantum motion of nanobeam resonances. 
Excitation of nanomechanical self-oscillations with nW 
absorbed power illustrates the sensitivity of the diamond 
nanobeams to small driving forces, and their nonlinear 
dynamical softening provides a glimpse of the changing 
stress within the nanobeam during large amplitude oscil¬ 
lations. This demonstration of optomechanical excitation 
of the diamond nanomechanical environment is a step to¬ 
wards optomechanical control of quantum electronic sys¬ 
tems such as nitrogen vacancies [TMn Ea |34]. It is 
promising for implementations of nanomechanical logic 
in diamond and can be extended to use the opti¬ 
cal gradient force for optomechanical nanobeam actua¬ 
tion, enabling excitation of higher frequency resonances 
of smaller structures Finally, the scalable nanofab¬ 
rication technique demonstrated here is widely applica¬ 
ble to diamond nanophotonic devices for sensing, non¬ 
linear optics, and quantum information processing, and 
can be easily adopted by researchers with access to stan¬ 
dard nanofabrication tools. Extending this approach to 


integrate an optical cavity into the waveguide optome¬ 
chanical system, for example by patterning the waveg¬ 
uide ends with mirrors [22], will allow studies of dissipa¬ 
tive optomechanical cooling of the nanobeam resonances 
in the unresolved-sideband regime m- 
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Appendix A: Fabrication process 

The fabrication process flow is outlined in Figure 2. 
A chemical-vapor-deposition-grown, (lOO)-oriented sin¬ 
gle crystal diamond optical-grade substrate (Element 
6) is cleaned in boiling piranha (3:1 H 2 S 04 :H 202 ) and 
coated with a 300-nm-thick layer of PECVD Si 3 N 4 as 
a hard mask. Next, 2 — 5nm of titanium is deposited 
as an anti-charging layer. Nanobeam structures with 
axes aligned along the (110) direction are patterned in 
ZEP520A electron-beam lithography resist and devel¬ 
oped at —15 °C in ZED-N50. The resulting ZEP pat¬ 
tern is transfered to the nitride hard mask using an 
inductively-coupled plasma reactive-ion etch (ICPRIE) 
process with C 4 F 8 /SF 6 chemistry. An anisotropic oxy¬ 
gen plasma ICPRIE step transfers the pattern to the di¬ 
amond, followed by a 160-nm-thick conformal coating of 
PECVD Si 3 N 4 to protect the vertical diamond sidewalls. 
A short anisotropic C 4 F 8 /SF 6 ICPRIE step clears the 
nitride from the bottom of the windows while keeping 
the side and top surface of the devices protected. To 
create suspended nanobeams, a quasi-isotropic oxygen 
plasma etch |SB| is performed at 2500 W ICP power, 0 
W RE, and elevated wafer temperature of 250 °C. This 
etch step is relatively slow, requiring ~ 5 hours to un¬ 
dercut the nanobeams studied here. The undercut etch 
rate was observed to increase at higher-ICP power, how¬ 
ever this option was not available for the devices fabri¬ 
cated for this report. An Oxford Plasmalab is used for 
all plasma etch steps. Finally, the titanium and nitride 
layers are removed by wet-etching in 49% HE, and the 
sample is cleaned a second time in boiling piranha. Note 
that adding an second vertical diamond etching step im¬ 
mediately prior to the quasi-isotropic etch is expected 
to reduce the necessary undercut time, as in the related 
SCREAM silicon process [43] . 

Appendix B: Evanescent coupling to diamond 
nanobeams 

Numerical prediction of the coupling coefficient k(/i) 
describing the interaction between the nanobeam and 



fiber taper waveguide modes can be obtained from the op¬ 
tical dispersion of the eigenmodes of the coupled waveg¬ 
uides (“supermodes”). This process is described below. 

The field propagating through the coupled waveg¬ 
uide system can be represented as a superposition of 
modes of the uncoupled optical fiber taper and diamond 
nanobeam. In the case of waveguides with two nearly 
phase-matched co-propagating (positive group velocity) 
modes, this evolution can be approximately described by 
coupled mode equations describing the field amplitude 
af_n(z) in the fiber and nanobeam waveguides, respec¬ 
tively, as a function of propagation distance z through 
the coupling region: 


daf 
dz 
da a 
dz 


-j{^i+ Kff)af- jKQn, 


-jiPn + Knn)an “ JKOf., 


(Bl) 

(B2) 


where /3f and /3„ are the A dependent propagation con¬ 
stants of the uncoupled fiber and nanobeam modes, re¬ 
spectively. Also included in this model are “self-term” 
corrections Kff,nn(/i) to /3f,n resulting from the modifica¬ 
tion of the local dielectric environment by the coupled 
waveguides. From Eqs. (IBII) and (B2), supermodes of 


the coupled waveguides, can be found with z dependence 
g-*/3±z^ where 








+ 


(B3) 


with /3f_n = /3f^n(A)-l-Kff_„n(ft,) [IS]. Supplementary Figure 
shows /3±(/i, A) calculated with a mode solver (Lumer- 
ical MODE Solutions). k(/i) and Avff.nn( fe) w ere then de¬ 
termined by fitting P±{h,X) with Eq. (B3), where nu¬ 


merically calculated /3f,n(A) of the uncoupled waveguides 
are known input parameters. The resulting values for 
k(/i) and Kff,nn(^)5 shown in Supplementary Fig. are 
found to decay exponentially with h. For simplicity, it 
was assumed that Kg = Knn during the fitting process. 
Note that while K{h) decays over a length-scale com¬ 
parable to the waveguide evanescent field decay, ng^nnih) 
are near-held terms which decay much more quickly, and 
do not signihcantly impact the predicted values of k(/i) 
for h > 100nm. For a given h, P±{X) can be converted to 
an effective index dispersion curve, n±(A) = /3±(A)A/27r, 
as shown in Fig. 1(c), which clearly exhibits an anti¬ 
crossing near /3f = /3n whose width scales with |k|. The 
excellent agreement between the numerically calculated 
P± and the semi-analytic model described by Eq. (B3) 


indicates that this model is suitable for predicting the 
waveguide coupling. 

The coupler response can also be pred icted from this 
coupled mode analysis. Solving Eqs. (Bl) and (B2| with 


boundary condition corresponding to unity input power 
to the hber taper (|af(0)p = l,|an(0)p = 0), results in 
the taper transmission T{X,h) given by Eq. (1) in the 
main text. 



FIG. 6. Propagation constants (/3±) of the even and odd 
supermodes of the coupled nanobeam and fiber taper waveg¬ 
uides, as a function of waveguide separation h, for varying A. 
Points indicated by open squares and circles where calculated 
using Lumerical MODE Solutions. Solid lines are fits using 
the coupled mode theory model described in this section. 



FIG. 7. Waveguide coupling coefficients as a function fiber 
taper and nanobeam waveguide separation (h). Each data 
point was obtained by fitting /3±(A) with the coupled mode 
theory model. Solid lines are single exponential fits to the 
data points. 


Appendix C: Measurement setup 

The optomechanical properties of the diamond 
nanobeam waveguides were studied by monitoring the 
optical transmission of a dimpled optical fiber taper po¬ 
sitioned in the near field of devices of interest. The dim¬ 
pled fiber taper was fabricated by modifying the proce¬ 
dure presented in Ref. |44j to use a cermanic mold for 
creating a dimple. Measurements performed in ambient 
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conditions used high-resolution (50 nm) DC stepper mo¬ 
tors to position the fiber taper. Vacuum (room temper¬ 
ature) and cryogenic measurements were performed in a 
closed cycle cryostat (Montana Systems Nanoscale Work¬ 
station) whose sample chamber is configured with stick- 
slip and piezo positioning stages (Attocube) for control¬ 
ling the sample and fiber taper positions. Before cooling, 
a turbo pump was used to evacuate the chamber to pres¬ 
sures in the 10“® Torr range. Room-temperature vac¬ 
uum measurements were performed in these conditions. 
At 5K, cryo-pumping reduced the chamber pressure to 
the 10“® Torr range. During the low temperature mea¬ 
surements, the fiber taper was positioned in contact with 
lithographically defined supports on the diamond chip to 
reduce coupling of vibration from the cryostat cooling 
stages to low frequency resonances of the optical fiber 
taper. These supports, visible in Supplementary Fig. 4, 
allow the fiber taper to be positioned in the nanobeam 
near-held without contacting the nanobeam. 

Two external cavity tunable diode lasers (New Focus 
Velocity 6700) were used to probe the hber taper trans¬ 
mission T over a wavelength range from 1475 - 1625 
nm. A New Focus 1811 photodetector (PDl) with a 
noise equivalent power =2.5 pW/-\/Hz was used to 
monitor the average (Pd) and huctuating {6Pd{t)) output 
power from the hber taper. A New Focus 1623 detector 
(PD2) was also used in some measurements of Pd- 

A Tektronix RSA5106A real time spectrum analyzer 
(RSA) allowed fast spectral analysis during the experi¬ 
ments, and recording of IQ time-series of the PDl output 
voltage V{t). All of the Sy{f) data presented here was 
generated in post-processing from V{t) data. By choos¬ 
ing a low center (demodulation) frequency (typically 0 
or 2.5 MHz), and sampling V{t) with a bandwidth ex¬ 
ceeding the nanobeam resonance frequencies (typically 5 
MHz) both ST{t) and T could be recorded by the RSA. 
To avoid damaging the RSA with a large DC input, a 
pair of bias-Ts (Minicircuits ZFBT-6GWB-I-) together 
with electrical attenuators were used to reduce the low 
frequency (< 100 kHz) components of V{t). The self¬ 
oscillation data in Fig. 5 was then acquired in a single 45 
s time-series while the hber height above the sample was 
stepped in 30 nm increments. Attenuation of the low- 
frequency signal was compensated for in post-processing. 


Appendix D: Optomechanical transduction 
sensitivity 


This section describes the procedure for predicting the 
theoretical measurement sensitivity of the optomechani¬ 
cal waveguide readout, and for extracting the experimen¬ 
tally observed measurement sensitivity from measured 
thermomechanical signals. Single-sided power spectral 
densities are used throughout. 


a. Theoretical sensitivity 

The displacement sensitivity of the waveguide- 
optomechanical system is determined by the minimum 
mechanical motion to actuate a signal larger than the 
noise hoor of the measurement apparatus. For direct 
photodetection of the optical power transmitted by the 
waveguide, the power spectral density, S'«*^(/), of the 
transduced signal from a mechanical displacement de¬ 
scribed by spectral density Sx{f) is 

si^\f) = SM) , (Di) 

where gti is the detector transimpedance gain, Pd is the 
detected power in absence of coupling (T = 1), and 
\dT/dx\ describes the optomechanical actuation of the 
coupler. Similarly, the measurement noise can be writ¬ 
ten as 


= + (D2) 


describes noise intrinsic to the detector, and is 
related to the detector’s noise equivalent power figure 
(^(det)) The contribution from pho¬ 

ton shot noise is given by 




2hwoTPd 


Vqe 


(D3) 


where rjqe is the detector quantum efficiency. Note that 
the impact of shot noise is affected by the operating 
point, T, of the coupler. This analysis does not con¬ 
sider optomechanical back-action, which is small for the 
optical powers used in the high sensitivity measurements 
presented here. 

To calculate the minimum detectable spectral density, 
S°, we require unity signal to noise ratio, 
resulting in 


SZ = 


2huJoTPd/riqe + S 


(det) 


(D4) 


5° has units of nm^/Hz and is related to the minimum 
detection sensitivity given in the text by s° = In 

the case of a co-directional evanescent coupler, dTjdx 
can be calculated from Eq. (1) in the main text. For 
perfect phase-matching (A/3 = 0), 


dT 

dx 


sin(2KL) 


dn dh 
dh dx' 


(D5) 


Alternately, dT/dx can be measured experimentally by 
recording T{h) and determining dT/dh. This approach 
was used together with Eq. ( |D4[ ) to generate the pre¬ 
dicted measurement sensitivity in Fig. 4, assuming that 
dh/dx = —1. 
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b. Thermomechanical calibration 


The observed thermal nanobeam motion was used to 
calibrate the measurement noise floor using a standard 
procedure described in, for example, Refs. [5^[Sn]- Ther¬ 
momechanical resonances in were fit using Sy{ijj) = 

5^"^ -I- GSl.^{u}) where is the single-sided power 

spectral density of a thermal oscillator. 


Sf{^) = 


AksTeCJn 


Qm ^ [(^ d” /Qm')^] 


(D6) 

and G is a constant determined by the transduction gain 
of th e measurement, as described theoretically by Eq. 
(D11, and treated as a fitting parameter for the purpose 
of the calibration procedure. Here ks is Boltzmann’s 
constant, is the operating temperature, and m is the 
effective mass of the resonance, as defined in Ref. [HT] . 
From the fit values, spectra can be converted from elec¬ 
trical (W/Hz) to displacement (m^/Hz) power spectral 
density: Sx = Sy/G. 



FIG. 8. Optical gradient force as a function of h, when 
the fiber taper is aligned along the nanobeam axis, and off¬ 
set 500 nm laterally, approximating the position in the self¬ 
oscillation measurements. 


Appendix E: Optical gradient force 


Dielectric objects in an evanescent held experience an 
optical gradient force. The optical gradient force between 
coupled waveguides can be predicted from n±{h) using 
the formalism of Povinelli et al. |62]- Following this for¬ 
malism, the force induced by power P± in each of the 
supermodes is given by 


F± 


P±Lc 


dh 


P±Lc 

c 



(El) 


Operating at phase matching with power Pi input into 
the hber taper, the power in the even and odd super¬ 
modes is P± = Pi/2, and the total optical gradient force 
on the nanobeam is 


„ _ PiLy / dn+ dn_ \ 
® “ 2c \dh ^ dh ) 


(E2) 


where positive (negative) Fg indicates a repulsive (attrac¬ 
tive) force. Supplementary Figure]^ shows the predicted 
Fg{h)/Lc Pi, for varying alignment of the hber taper with 
the center axis of the nanobeam. Corrections due to the 
curvature of the dimpled hber taper are not considered 
explicitly, and are assumed to be accounted for by the ef¬ 
fective coupler length Ly extracted from the experimental 
measurements. Given this approximation, Fg is constant 
over the interaction length of the ideal two-port waveg¬ 
uide coupler, and vanishes outside of the coupling region. 


Appendix F: Nanobeam photothermal and nonlinear 
dynamics 

Waveguide-optomechanical coupling was observed to 
induce nanobeam self-oscillations, as shown in Fig. 5 of 


the main text. Here we analyze this effect by modeling 
the nanobeam as a nonlinear harmonic oscillator inter¬ 
acting with the optical held of the coupled waveguides 
through a dynamic photothermal force and an instanta¬ 
neous optical gradient force [n i3il[5?[|M] . 


c. Nanobeam equation of motion 


The dynamics of the nanobeam resonance driven by a 
stochastic thermal force, Fs{t), and coupled to the optical 
hber taper, are approximately described by the following 
equation of motion: 


X + + 0 : 32 :^ + C( 2 X)X = 

— {Fg{x) + Fptix,t) + Fs{t)) , (FI) 

m 

where x{t) is the displacement amplitude of the 
nanobeam mechanical resonance of interest, dehned rela¬ 
tive to the position of the undriven nanobeam. Changes 
in X modulate the nanobeam and the hber taper spacing 
h, and for the case of the Vi resonance, h{x) = hg — x{t) 
if the hber taper is positioned ho above the center of the 
nanobeam. Optomechanical coupling arises from the de¬ 
pendence of the optical gradient force, Fg{h{x)), and the 
photothermal force, Fpt{h{x)), on nanobeam position. In 
addition, uim varies due to internal strain resulting from 
optically induced changes in static dehection, x, and lo¬ 
cal nanobeam temperature relative to the environment, 
0. For sufficiently large |x(t)|, the nanobeam response 
becomes nonlinear, as described by 02 and a^. 
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d. Thermal effects 


Local heating of the nanobeam by waveguide opti¬ 
cal absorption occurs on a timescale determined by the 
nanobeam geometry and material properties. The ther¬ 
mal dynamics of the nanobeam, defined by the maximum 
temperature change, 0, relative to the operating temper¬ 
ature, are assumed to follow 

^ = -KtS + RL.CPn, (F2) 

at 

where r = I/k* is the nanobeam thermal time-constant 
and Pn is the optical power coupled into the nanobeam. 
In absence of delayed optical feedback (i.e., a cavity), P„ 
instantaneously follows x, and is given by 


= (1 - T{h{x)))P„ (F3) 


with T{h) described by Eq. (1) in the main text. Here 
C is the per-unit length absorption coefficient of the 
nanobeam waveguide, and R is the heating power per 
unit of absorbed optical power. Ideally, C, is determined 
by the material properties of diamond, but in nanopho¬ 
tonic devices it can be modified by imperfect surfaces. 
The effective nanobeam waveguide optical interaction 


length is given by Li{h) = \an{z)\^dz/\af (Q)\^, and 
represents the distance over which light propagates in 
the nanobeam. It can be calculated from Eqs. (Bl) and 

Local heating of nanobeam waveguides induces deflec¬ 
tions in the nanobeam position. We write the corre¬ 
sponding photothermal force as 


Fpt{x,t) = Pe{x,t)'^, (F4) 

where P is the force per unit absorbed power (Lit^Pn) in 
steady state. The resulting change in nanobeam deflec¬ 
tion is Sxpt = Fpt/k where k = ui'^m is the nanobeam 
spring constant. 

Local heating and accompanying thermal expansion of 
the device also modifies ujm- This effect is represented by 


doJm _ ^ 

de ~^*R 


(F5) 


where Ct is a constant related to the elastic properties of 
the nanobeam, and has units rads“^W“^. 


e. Static response of compressed nanobeams 

In general, Ct and P sensitively depend on both the 
nanobeam elastic properties and geometry, and the inter¬ 
nal residual stress acting on the device. The nanobeams 
used in the self-oscillation studies have significant inter¬ 
nal stress, manifesting in smaller Wm than expected from 
their nominal dimensions, and can be in a buckled ge¬ 
ometry as shown in the SEM image in Supplementary 


Fig.0 To model nanobeam behavior in the presence of 
compressive stress and buckling, together with imperfect 
nanobeam shape and clamping points, we consider both 
an approximate analytic model, and finite elements sim¬ 
ulations. 

The behavior of an ideal beam under axial compressive 
loads has been widely analyzed [5S] . Eor compressive ax¬ 
ial load Fi, uJm of a nanobeam with maximum deflection 
X can be expressed as 




F,, 




= «)^ l-^+Sx 


AE 


Pc (2F)2 


(F6) 


where is the resonance frequency of an ideal un¬ 
loaded nanobeam, A is the nanobeam cross-sectional 
area, E is Young’s modulus and F^ is the critical buckling 
load. Equation (F61 is equivalent to the model presented 


in Ref. [55], and reveals the interplay between deflec¬ 
tion, axial load, and stiffness. For example, we see that 
dujm/dx oc X, which is a result of more efficient conver¬ 
sion of transverse actuation to axial strain with increas¬ 
ing nanobeam curvature. Locally heating the nanobeam 
modifies Ft —>■ Ft + rjeEQ, where e is the thermal ex¬ 
pansion coefficient of the nanobeam and 0 < ry < 1 is a 
geometric factor which accounts for non-uniform heat¬ 
ing distribution of the nanobeam. In general, F) and 
X are not independent. However, in an ideally buckled 
nanobeam 1651. we can show that 


dx P Kt _ ijewdE 

dO k R 2muj^L’ 


(F7) 


indicating that photothermal deflection can be enhanced 
in buckled nanobeams with \x\ > 0. 

Although the ideal nanobeam buckling model is in¬ 
structive, it fails to reproduce experimental features such 
as deflection for axial load below the critical buckling 
load, and effects related to imperfect elastic clamping 
points and beam deformation |5S]. To predict P and 
Ct while including nanobeam non-idealities, we used fi¬ 
nite element ANSYS software to simulate ojm and x as 
a function of axial and transverse loads, and for speci¬ 
fied absorbed power. The results, summarized in Sup¬ 
plementary Fig. 10 and Table [l] indicate that significant 
enhancement of P compared to an unloaded nanobeam 
are expected. 

The simulations were conducted as follows. The simu¬ 
lated structure consists of a nanobeam with dimensions 
Lxwxd = 80 X 0.48 x 0.25/im^, and includes the 
surrounding diamond chip. Clamping point geometry 
was found to significantly affect the simulated nanobeam 
properties. In fabricated nanobeams, the undercut pro¬ 
cess results in relatively complex clamping point geome¬ 
try. Here the clamping points were modeled with a trian¬ 
gular vertical profile roughly approximating that of fab¬ 
ricated structures, extending 0.5 /rm on the nanobeam 
bottom surface and 1.0/rm vertically along the under¬ 
cut sidewall, as shown in Supplementary Fig. |10[ a). uJm 
and X were then calculated as a function of Fi, as shown 
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FIG. 9. SEM image of typical nanobeams observed to self-oscillate. Dashed straight lines are guides to emphasize the direction 
of nanobeam buckling. 


in Supplementary Fig. |10[ b). Agreement between sim¬ 
ulated and experimental cUm was realized at two values 
of Fi, corresponding to pre- and post-buckled nanobeam 
states, with x = —9 nm and —123 nm, respectively, where 
negative x indicates buckling down. 

F and Ct were estimated by simulating changes to x 
and uJm when power Pabs is uniformly absorbed across 
half of the nanobeam. During these simulations the bot¬ 
tom surface of the diamond substrate was fixed at con¬ 
stant temperature, and only conductive heat loss was 
considered. The corresponding temperature distribution 
was used to predict i? and kj, as summarized in Table 
|lj The photothermal results are shown in Supplementary 
Fig. [loj^c), which plots the changes Ax and Acum in de¬ 
flection and frequency, respectively, for Pabs = 1 /^W. 
These results clearly illustrate the sensitivity of pho¬ 
tothermal effects on compressive stress and deflection, 
and indicate that F varies by over two orders of magni¬ 
tude depending on the compressive axial loading of the 
nanobeam. 

Simulations of AcUm when a vertical transverse load 
mimicking the optical gradient force is distributed along 
the coupling region of the nanobeam (length Pc) were 
also performed, and are shown in Supplementary Fig. 
[I0|(d). These results show that doJm/dFg varies by nearly 
three orders of magnitude depending on the compressive 
axial load. 


/. Dynamics: small-amplitude oscillations 


For small mechanical oscillation amplitude, the power 
in the nanobeam waveguide can be approximated by 


Pn = PniK) - Px 


dT I 
dx \ ho 


{x{t) — x). Inserting this into 


Eq. (F21 and only retaining terms that are linear in x 
allows the Laplace transform of Eqs. (FI I and (F21 to be 
combined into a single linear equation: 


2 / \ .^ni^ / \ I 2 / \ 

— OJ X(UJ) — 1-— - X[UJ) -I- UJ^X{UJ) = 

Fs{u}) fdFg dFpt wT-bl \ x{u}) 
m \ dx dx + 1 / ^ 


where constant force terms resulting in changes to static 
nanobeam deflection x have been left out for clarity 
but are included implicitly in Rearranging 

terms reveals the optomechanical renormalization of the 
nanobeam dynamics. 


- U}‘^x{uj) - i-f'^ujx{uj) + uj'^xiuj) 


Fsjoj) 

m 


with 

^ ^ _ 1 dFptjhp) 1 _ dFg{ho) 1 

1 -I- dx k dx k 

7m _ 1 I ^ dFpijho) 

7m ™ 1 -I- dx k ’ 


(F9) 


(FIO) 

(Fll) 
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FIG. 10. (a) Cross-section of nanobeam used in ANSYS 

finite element simulations, showing the “arched” clamping 
point geometry, (b) Finite element simulations of resonance 
frequency and nanobeam deflection as a function of axial 
compressive stress. Dashed lines indicate value of Fi where 
simulated cdm approximately matches the experimentally ob¬ 
served value, (c) Simulated change in resonance frequency 
and nanobeam deflection for Pabs = 1 of absorbed opti¬ 
cal power, (d) Simulated change in resonance frequency from 
a Pg = 1 pN transverse load applied across the coupling re¬ 
gion of the nanobeam (length L^)- Nanobeam dimensions are 
Lxwxd = 80 X 0.48 x 0.25/rm^, as in the self-oscillation 
measurements. 


where dFpt/dx = —Fd^LiPidT/dx. Static effects modify 
the unperturbed resonance frequency according to 


= < + CtCL, Pn{h) + ^Fg{h). (F12) 

dtg 


Equations (FIO), (Fill and (F12) illustrate that combi¬ 


nations of mechanical softening or hardening, and ampli¬ 
fication or damping are possible. Higher order thermoe¬ 
lastic effects [S3] are predicted to be small and are not 
included in this analysis. 


grows, and nonlinear contributions to the system dynam¬ 
ics become significant. These nonlinear modifications 
originate from mechanisms intrinsic to the nanomechani¬ 
cal device geometry, or to the optomechanical response of 
the system, and are characterized here by non-zero 02,3 
and d^T/dx'^ (n > 1), respectively. Two nonlinear fea¬ 
tures observed in Fig. 5 are a softening in uj^n at the onset 
of self-oscillation, and frequency harmonics in the self¬ 
oscillation region. While the later effect is significantly 
affected by the nonlinear response of the optomechani¬ 
cal system, the analysis of Zaitsev et al. [M] shows that 
the optomechanical nonlinearity plays a negligible role 
in softening Um for the optomechanical system studied 
here. Rather, this softening is dominantly due to geo¬ 
metric nonlinearities of the deflected nanobeam. 


Nonlinear coefficients 02 and 03 in Eq. (FI), can be de¬ 


rived from the Euler-Bernoulli equation for a nanobeam 
with static deflection x^ as in Refs. [67llM]- Using an ap¬ 
proximate ansatz for the static nanobeam shape, expres¬ 
sions for Q!2 and 03 can be derived |5S]. For oscillations 
about x^ if we group the nanobeam deformations into 
time dependent and time independent parts according to 

t) = {x(t) +x)^{l — cos{2'xl/L )), (F13) 

where I is the the coordinate running the length of the un¬ 
deflected nanobeam, the Euler-Bernoulli equation yields 
Eq. (FI I, with nonlinear coefficients 


E , 

f2Tr\^ 3 

(F14) 


\T) 8’ 

E 1 

f27ry 3 

(F15) 

^ 1^ ( 

vTJ 8’ 


where p is the density of the nanobeam material. Note 
that while 02 vanishes in a straight nanobeam, it is non¬ 
zero in deflected nanobeams (|s| > 0). _ 

Assuming that the solution to Eq. ( |F1| ) may be writ¬ 
ten as a combination of harmonic functions, the method 
of successive approximations shows that to first 

order the fundamental frequency of oscillation is: 


' = UJm{x) + 


Ulm{x) 



5 

12w^(x) 



(F16) 

(F17) 


where v is the amplitude of oscillation and a is the effec¬ 
tive nonlinear frequency shift coefficient. The first term 
in brackets is the well known Duffing frequency modifica¬ 
tion, whereas the second term results from nonlinearities 
induced by static deflection of the nanobeam. 


g. Nonlinear dynamics: large-amplitude oscillations 

When approaches zero, the amplitude of mechan¬ 
ical oscillations about the static (i.e., buckled) position 


h. Parameter estimation and comparison with experiment 

To compare the model described above with experi¬ 
mentally measured as shown in Fig. 5(c), a combi¬ 

nation of known, simulated and fit parameters were input 
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FIG. 11. Predicted generated using the model in Eq. 

(FlOl and the parameters in Table 


into Eq. (Eli), as summarized in Tableland described 
below. 

For a given coupler operating condition, the photother- 
mal force, and the resulting 7'(/i) described by Eq. (Fill, 
scales linearly with Neither or ^ are known a 

priori. However, using finite element simulations to de¬ 
termine F, as well as other parameters such as Kt and 
i?, allows C to become the sole fitting parameter when 
comparing experimental and predicted values of ^'{h). 

To determine F from finite element simulations, it is 
necessary to determine whether the nanobeam is in a 
pre- or post-buckled state. Insight into the buckling 
configuration of the nanobeam is provided by the ob¬ 
served nonlinear softening at the onset of self oscilla¬ 
tions, which is directly proportional to (Fig. 5(d)). 
From this data and Eq. (F17), the nonlinear coefficient 
a given in Tabl e [I| c an be measured. Together with 
Eqs. (F14) and (F15), corresponding nanobeam deflec¬ 
tion amplitude of |3;|= 98nm is inferred, in good agree¬ 
ment with X = — 122nm predicted from finite element 
simulations of the nanobeam in its post-buckled config¬ 
uration {FijA ~ 37MPa), as shown in Supplementary 
Fig.[^b). When the nanobeam is in this post-buckled 
configuration, finite element simulations shown in Sup¬ 
plementary Fig. [lO|(c) indicate that F = —26pN//xW 
(Fjk = —0.36 nrn7/rW). Negative F and x indicate that 
the photothermal force and nanobeam deflection, respec¬ 
tively, are in the down direction. 

Inputing F and other finite element simulated pa¬ 
rameters summarized in Table [l] into Eq. (Fll), good 
agreement between predicted and observed q^h) was 
found for ^ ^ 0.12cm“^. The corresponding optical ab¬ 
sorption rate can be described by quality factor Qo = 
2'nnglC,\ ~ 6.6 x 10^, where the group index rig ~ 2.0 
of the nanobeam is predicted from numerical simulation. 
This absorption rate is smaller than combined absorption 
and radiation loss rates in other single crystal diamond 
nanophotonic structures [i^ . 

Given the value of C obt ained from fitting 7'(/i), 
predicted from Eq. ( |F10 ) can be compared with mea¬ 
surements. This is shown in Supplementary Fig. [m in 
which the predicted was generated with duJm/dFg 

as a fitting parameter, and other parameters set as in 
the model for 7'(/i) and listed in Table |T] Static and dy- 


Parameter 

Input value 

Units 

Source 


1.4 


FEM simulation 

R 

0.11 

K/pWfis 

FEM simulation 


680 

kHz 

FEM simulation 

m 

10 

Pg 

FEM simulation 

X 

-122 

nm 

FEM simulation 

0.1211 

-16 (-28) 

Hz/nm^ 

Eit (FEM) 

Ct 

1.0 

kHz/^W 

FEM simulation 

C 

0.12 

cm“^ 

Fit 

F/k 

-0.36 

nm/p,W 

FEM simulation 

F 

-26 

pN^W 

FEM simulation 

Fo 

Sup. Fig.| 


Optical mode solver 

dujf-n 1 dFg 

-7163 (-377) 

rad/s pN“^ 

Eit (FEM) 

R 

300 

pW 

Experimental parameter 

Qm 

25000 


Experimental parameter 


TABLE I. Parameters input to model for 'y'mih) and 
Values in brackets are predictions from simulations, and are 
included for comparison with values determined from experi¬ 
mental fits. 


namic thermal effects are found to be significantly smaller 
than the maximum experimentally observed shift to LOm- 
Rather, the monotonic decrease in a;(„, which becomes 
significant for h < 200 nm, follows an h dependence con¬ 
sistent with static tuning by the attractive optical gradi¬ 
ent force Fg(h) pulling up on a down-buckled nanobeam. 
This is in contrast to the static thermal tuning described 
by Ct, which is proportional to the power coupled into 
the nanobeam, and is expected to decrease in magnitude 
with decreasing h for h < 200 nm. 

There is good agreement between the model and exper¬ 
imental observation in Supplementary Fig. 0 however 
the fit value for dWm/dFg is larger than expected from fi¬ 
nite element simulations (Supplementary Fig. |10[ d)) and 
predictions of Fg{h) (Supplementary Fig. 1^. As shown 
in Supplementary Fig. , dujm/dFg was found in sim¬ 
ulations to be highly variable, spanning over three orders 
of magnitude, depending on compressive stress and re¬ 
sulting buckling configuration. It is possible that im¬ 
perfect nanobeam and clamping point shape result in 
an enhanced sensitivity [Si Eg. For example, including 
the triangular “arched” clamping points to approximate 
the fabricated structure, as shown in Supplementary Fig. 
[^a), enhanced duim/dFg by a factor of ^ 2 compared 
to the case of ideal clamping points. Other effects not 
included here include reflection of light from the end of 
the waveguide resulting in enhance optical interactions, 
dynamic and static displacement of the fiber taper due 
to optical forces EH, low frequency fiber vibrations and 
possible parametric driving for small h, breakdown of 
the coupler two mode representation and the influence of 
higher-order modes for small h, and short range effects 
such as the Casmir force m- 

In conclusion, the analysis presented here serves to il¬ 
lustrate the influence of compressive loading and buck¬ 
ling on the photothermal response of nanobeams, and to 
show that significant photothermal forces are present at 
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relatively low optical absorption levels. Note that imper¬ 
fect nanobeams have significantly altered pre- and post- 


buckling behavior when subject to an compressive load 
m, and taking into account these non-idealities is neces¬ 
sary to more accurately predict the nanobeam behavior. 
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